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Abstract 

It is possible to fabricate mesoscopic structures where at least one 
of the dimensions is of the order of de Broglie wavelength for cold elec- 
trons. By using semiconductors, composed of more than one material 
combined with a metal slip-gate, two-dimensional quantum tubes may 
be built. We present a method for predicting the transmission of low- 
temperature electrons in such a tube. This problem is mathematically 
related to the transmission of acoustic or electromagnetic waves in a two- 
dimensional duct. The tube is asymptotically straight with a constant 
cross-section. Propagation properties for complicated tubes can be syn- 
thesised from corresponding results for more simple tubes by the so-called 
Building Block Method. Conformal mapping techniques are then applied 
to transform the simple tube with curvature and varying cross-section 
to a straight, constant cross-section, tube with variable refractive index. 
Stable formulations for the scattering operators in terms of ordinary dif- 
ferential equations are formulated by wave splitting using an invariant 
imbedding technique. The mathematical framework is also generalised to 
handle tubes with edges, which are of large technical interest. The numer- 
ical method consists of using a standard MATLAB ordinary differential 
equation solver for the truncated reflection and transmission matrices in a 
Fourier sine basis. It is proved that the numerical scheme converges with 
increasing truncation. 



1 Introduction 

In the search for faster computers critical parts are becoming smaller. Today, 
it is possible to build mesoscopic structures where some dimensions arc of the 
order of the de Broglie wavelength for cold electrons. Often the electron motion 
is confined to two dimensions. Consequently, it may be necessary, at least for 
some computer parts, to include quantum effects in the design process. 

A large number of studies, devoted to such quantum effects, have been car- 
ried out in recent years and a review is given by Londegan et al R]. Many 
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investigations aim at understanding the physical properties of a particular quan- 
tum tube rather than developing reliable mathematical and numerical methods 
that can be used in a more general context. The research has given valuable 
knowledge on the physical behaviour but also reports on the limitations of the 
methods used. For instance, Lin & Jaffe |^ report that a straightforward match- 
ing at the boundary of a circular bend does not converge, demonstrating the 
numerical problems with such a method. An illposedness is present in quantum 
tube scattering and some type of regularisation is therefore required to avoid 
large errors. Often, the tubes have sharp corners to facilitate manufacturing 
but also to enhance quantum effects. The presence of corners with attached 
singularities requires special treatment. 

Scattering of electrons in quantum tubes, see figure 1, is theorywise related to 
the scattering of acoustic and electromagnetic waves in ducts. Nilsson treats 
a general method for the acoustic transmission in curved ducts with varying 
cross-sections. Wellposedness, i.e. stability, is achieved in an asymptotic sense. 
The mathematical framework guarantees consistent results and allows for sharp 
corners and a proof for numerical convergence is given. We set out to present a 
quantum version of the results of Nilsson Q . In this way the problems reported 
on convergence and on inconsistent mathematical results would be resolved. 



Figure 1: Two-dimensional quantum tube 



The paper is organised as follows. An introduction to scattering in quantum 
tubes is given in section 2 and a mathematical model is formulated in section 3. 
The Building block Method which is a systematic method to analyse complicated 
tubes in terms of results for simple tubes is also briefly described. Then in 
section 4 the scattering problem for the curved tube with varying cross-section 
and constant potential is reformulated to a scattering problem for a straight 
tube with a varying refractive index. The solution to this problem is presented 
in section 5 and a discussion on numerical methods are also given. 
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2 Tubes in quantum heterostructures 



A schematic view of a quantum hctcrostructure is shown in figure 2 following Wu 
et al. Q Electrons are emitted from the n-type doped AIGaAs layer, migrate 
into the GaAs layer and stay close to the boundary to the AIGaAs layer. In 
this way a very narrow layer of electrons which are free to move in a plane is 
formed. Nearly all the electrons in this two-dimensional gas are in the same 
quantum state. By applying a negative potential on the metal electrodes on the 
top of the heterostructure in figure 1, the electrons are banished from the region 
below the electrodes. For relatively low voltages, the effective potential in the 
tube for one electron is close to the square-well potential. As a consequence 
the electrons in the two-dimensional gas are further restricted to a tube that in 
form is a mirror picture of the gap between the two electrodes. This quantum 
tube links the electrons between the two two-dimensional gases on both sides of 
the strip formed by the electrodes. 




Figure 2: Schematic picture of heterostructure and split-gate structure. 



3 Mathematical model 

Consider a two-dimensional tube with interior f2' according to figure 1. The 
boundary F' consists of two continuous curves, F^ and F'_, which are piecewise 
C^. The upper boundary F^ can be continuously deformed to F'_ within fl'. 
Outside a bounded region the duct is straight with constant widths a and b, 
respectively. These terminating ducts are called the left and the right termi- 
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nating duct or L and R for short. We use stationary scattering theory for one 
electron in an effective potential, with time dependence exp(—iEt/h), assuming 
that the wave function ip satisfies the time-independent Schrodinger equation 
Alp + k'^ip = in r2',where fc^ = 2m*E/h and m* is the effective mass Q. 
Usually is called energy. The effective potential is assumed to be a square 
well meaning that ip\-p, — 0. 

In a tube with constant cross-section the harmonic wavefunction ip can be 
uniquely decomposed in leftgoing and rightgoing parts hy ip — + xp^ . Super 
indices " -|- " and " — " indicate rightgoing or plus and leftgoing or minus waves 
respectively. Let and ip~^ be known incoming waves in the terminating 
ducts. V'j'ti is present in the left and ip'^ in the right one. Let us write 

V' = i^tn + R^i^tn + r-V-r™ in L , . 

^ = + i?-Vr„ + in R ' ^ ^ 

where for example the last two terms in (|^a) are minus waves and the equa- 
tion defines the left reflection mapping i?+ that maps the incoming wave to 
an outgoing one in L. The scattering problem consists of finding the mappings 
i?+,T~,i?~ and T+ as functions of energy for a given duct. In summary we 
have 

AV- + fc^V- = in f]' 

V'+ = ^tn in L ■ ^ ^ 

V'" = V-j^ in R 

There is always a solution to (||), and except for a discrete number of eigenen- 
ergies fc^ = — 1,2,3,..., the solution is unique. |^ When fc^ — kj, an 
eigenenergy, there exists a solution without incoming but with outgoing waves. 

The use of the Building Block Method or transfer matrix formalism 
is very efficient for the solution of scattering problems. In this method a tube 
with a complicated geometry is divided into two parts usually where the tube 
is straight. These two parts are converted to the type shown in figure 1 by 
extending the terminating tubes to infinity. A sub tube for the tube shown in 
figure 1 originates from the left part and is depicted in figure 3. The Building 
Block Method gives a procedure for calculating the mappings R^, , , and 
T"*" for the entire tube in terms of the corresponding scattering properties for 
the sub tubes. This procedure can be repeated to get several sub tubes. Rather 
than using a general numerical package for conformal mappings we have for 
the calculations in this paper employed the Schwarz-Christoffel mapping for a 
duct with corners and rounding the corners using the methods of Hcnrici 
Required analytic integrations are performed in MATHEMATICA. 

We recall the standard duct theory in a form that illustrates the illposed- 



4 




0.4 
0.2 

-1 -0.5 0.5 1 

Figure 3: Sub- tube with interior O' and upper boundary F^and lower boundary 
F'_. 6^ = 0.6. 



ness of the problem and we have 

oo oo 

V, = V+ + V" = E ^;te-">„(y) + E A;;e--"Vn(2/), (3) 



n=l 



with ipniu) = sin(nny / a) and a„ = ^k"^ — n^TT^/a^, Im a„ > 0. It is convenient 
to define the operator Bq by 

Bof = J2'^=l OinfnVn, /^-j 
fiy) = J2n=l Otnfn^n{y) 

We find that 3^=8^ + and dx^J^ = ii-BoV-'^- The initial value problem, 

V'+(0) = Vo, ^' 

is illposed for a; < 0, but not for x > 0. If an attenuated plus wave is marched 
to the left an exponential growth is found. To avoid the illposedness, tp is 
decomposed and the plus waves arc calculated by marching to the right and 
minus waves in the opposite direction. 



4 Reformulated scattering problem 

To be able to use powerful spectral methods it is advantageous to transform 
the tube to a flat boundary. It is enough, according to the Building Block 
Method, to consider the scattering in the sub tubes and wc restrict ourselves to 
the first part as shown in figure 3. One way of transforming the tube is to use a 
conformal mapping w{() transforming the interior Cl' of the tube with variable 
cross-section in the ( = x + iy plane (figure 3) to the interior Q of a straight 
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Figure 4: ij,{u,v) in the straight duct. Parameters as in figure 3. ^ ^ is the 
refractive index. 



tube with constant cross-section in the w — u + iv plane. The straight tube is 
described by -co < u < oo, < w < a. 
Introducing (j){u, v) — "tp^x, y) we get 

i>{u,0) = <j){u,a)^0,ueR ' ^ ' 

with B^{u) = + k'^fj.{u,v) and /i = \d(/dw\'^. ii{u,v)~^ can be denoted as a 
refr'active index for the straight tube. In figure 4, ji related to the simple tube in 
figure 3 is depicted. The factor v) is asymptotically constant at both ends 
of the tube or more precisely ~ /i±-(-0(e^l^"l), u — > ±oo with /i_ = 1 

and ^+ = {b/aY- 

We use a first order description and rewrite (^) as 



duc^ ) \ -B' J \ du4> 

To avoid illposedness the decomposition (j) = 0+ + 0^ is introduced which 
must be identical to the corresponding decomposition (|^) in regions where /i 
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is a constant. The new state variables {(f)'^,(f) ) are introduced via the linear 
relation 

( ic ) ( r ) ■ 

Solving (H) for 0+and (f>^ and taking the w-derivative and using (0) we find that 



where 



a = i[(a„C-i)C + iC^^B^ + iC] 



/3=i 
7-1 



-(a„c-i)c - ic-ij52 + ic 

S =Y[(a„C-i)C - iC-^B^ - iC]' 



(10) 



To generalize the concept of transmission operators we make them u-dependent, 
using a similar notation as Fishman [nOl : 



(11) 



assuming that ui < U2. and suppressing the explicit w-dependence. It is assumed 
for ( pT| ) that the scattering problem has a unique solution or that homogenous 
solutions are removed. A homogenous solution is usually called a bound state. 

Next we find a differential equation for the scattering operators r+(u2, ui), 
i?~(ui,U2), R'^{u2,ui) and T~(ui,U2) in ( pi] ) using the invariant imbedding 
technique ]lO| . It is required that the incoming wave from the right, (/)^(u2), 
is vanishing. Then put ui = u, find du4>~{u) from (pT|), use ( pi] ) once more to 
obtain 

duR^{u2,u) = 7 + SR+{u2, u) - i?+(u2, u)a - R+ {u2,u)(3R+ {u2,u), (12) 

In a similar manner we get 

duT^iu2,u) = -r+(w2, u)a - T^{u2,u)PR'^{u2,u). (13) 

The stability properties of ([l^ ) and ( p^ ) are of central importance. In the flat 
regions where B = B^ or _B_ we have C = B and duC^^ — implying that 
/3 = 7 = and a = -J = iS. Similarly (|l|) and (|l|) reduce to a„X+ = 
—iBX^, — i?+ or T+, equations which are well-posed for marching to the 
left. The initial values to accompany ( p^ ) and ( p^ are i?+(it2,M2) = and 
T^{u2,U2) ~ I, where / is the identity operator. 

We choose C ~ i?_ + f{u){B+ — i?_) that is independent of v. Here / is 
increasing and smooth with lim„^_oo/(it) — 0, and lim„^oo/(w) — 1- 
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5 Solution of the scattering problem 



For the numerical solution of the scattering operator we expand in a Fourier 
sine series and /x in a Fourier cosine series: 



(14) 



where ^n(w) = cos(n7r/a). Using the notation (j) — (00, 01, ...)'^ we find that 

f^+BHu)^{u)^0. (15) 
The matrix elements of B^(u) are given by 

B^(u)„m TrhMm+n(") - fhn-n{u) - fim + Ai«-m(u)] ^Sn,n, n, TO = 0, 1,2, . 

(16) 



and it is understood in (|16|) that /ii(u) = for negative /. 

For the tube in the physical ^— plane we require that locally both the poten- 
tial and the kinetic part of the energy are finite, that is both dxdy < oo 
and \Vtjj\'^ dxdy < oo for all finite regions X inside the tube. We say that 
tjj belongs to the Sobolev space Hj^^ meaning that tp and its first derivatives 
are locally square integrable. Transformed to the straight duct the local finite 
energy requirement means J^j \4>\^ ^dudv < oo and Jjj \ W(f>\^ dudv < oo for all 
finite regions U inside the tube. For a smooth boundary is more regular, and 
also the second derivatives of (j) are square integrable, that is G ^loc- I* follows 
from the theory of Grisvard ||l^ that also the second derivatives of 4> are square 
integrable, which means that G H^^^. According to a graph theorem |Q G 
^loc implies that 0(u, •) € H^/^(0,a), meaning that up to 3/2 derivatives are 
square integrable. To interpret this regularity with fractional derivatives we 
define, following Taylor jl^, the function space 

D, = lfeL^{0,a):f2\fn\'{l + ny<^\,s>0, (17) 

I n=0 J 

with / = J2'^=i fn^n and /„ = (/, ipn)/ {ipm Vn)- IS a Hilbert space with the 
norm 



(,/,/) = El/"l'(l + (18) 



Taylor |T| shows that Dq =L2(0,a), Di =Hi^(0,a), D2 =H2(0, a)nHj^(0, a) and 
that 9uDs = Ds_i, s > 1. In this terminology we have that for a smooth 
boundary 0(u, •) e D3/2. 
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The operator is self-adjoint on 03/2- Thus, we may define B± by 



B±f = Y, v/fc2Ai± - n%Va2/«<Pn, (19) 

assuming that the branch Im > of the square root is taken. It is clear that 
T+, ii+ and T~ are mappings D3/2 ^03/2 and B±: Dg ^.s-i, s>l. 

For tubes with edges in the C— duct things are a little more complicated. 
With no restriction on the sharpness of the edges we cannot improve that 
4> S implying 4>{u,-) eDi/2- Then, as an intermediate step in our cal- 

culations B±(j) should be in the space D_i/2- Such a derivative must of course 
be interpreted as a distribution. However, the end result, i.e. scattered wave 
function belongs to D 1/2- To generalise we define by duality for positive s 

D_, = 1-9;^ f{v)g{v)dv < 00 for all / e D,| . 

Multiplication byy^ is an operator D1/2 — *■ D-1/2 and if s > 1/2 we have 
the following mapping properties: B± : Dg — s- Ds-i,dv : Dg — > Ds„i, and T+, 
, i?+ and are mappings -^D^. 

The equations ( p"^ - p^ ) can only in very special cases be solved in a closed 
form. Therefore some type of numerical scheme is used. Generally a numerical 
method cannot give uniform convergence for the entire space Dg. In a practical 
application it is usually sufficient to know the effect of the scattering matrices 
on the lowest eigenfunctions, the first A'o say. A practical method is therefore 
to truncate the matrix representation of ( |l2|) - ( [l3| ) to >> A'o and solve 
the finite-dimensional ordinary differential equation with a standard numerical 
routine. Nilsson |3| proves that such a procedure converges when N ^ 00. 

Presently, numerical results are not available for the quantum tube scat- 
tering. However, Nilsson presents results for the acoustic case where the 
Neumann rather than the Dirichlet boundary condition applies. He reports 
that for the lowest order reflection coefficient = 1, i.e. a scalar solution, is ac- 
curate up to fca = 1.5, Af = 2 gives a good and N — 5 gives a perfect discription 
up to fca = 6. Energy conservation holds for all N. 
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